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Abstract 

' — ^  A  derivation  of  surface-dressed  optical  Bloch  equations,  involving 
a  treatment  of  surface-reflected  photons  and  a  surface  plasmon 
resonance,  is  presented  for  a  collision-damped  two-level  atom  near 
or  adsorbed  on  a  metal  surface.  Effects  of  the  laser  bandwidth  are 
Included  by  means  of  a  phase-diffusion  model  for  the  driving  field. 

In  the  weak-fleld  or  large  detuning  limit,  the  population  Inversion 
and  resonance  fluorescence  spectrum  are  obtained  analytically.  These 
quantities  along  with  the  surface-induced  phase-decay  constant  of  the 
adatom  show  strong  oscillatory  behavior  as  a  function  of  the  adatom- 
surface  distance.  <—■ 


♦Permanent  address:  Institute  of  Physics,  Academia  Slnica, 
Beijing  (Peking),  The  People's  Republic  of  China. 


I.  Introduction 

Recently  there  has  been  considerable  Interest  In  the  Interaction  of 
electromagnetic  radiation  with  surfaces.  The  experiments  which  have  moti¬ 
vated  this  present  paper  are  those  demonstrating  the  high  Intensity  of  light 
Inelastlcially  scattered  from  molecules  adsorbed  on  metal  surfaces*  and  the 

surface-enhanced  luminescence  of  an  excited  molecule  fluorescing  near  or  on 

2-4 

metals  and  dielectric  fibers. 

To  Investigate  the  basic  behavior  of  the  surface-enhanced  resonance 
excitation  of  a  laser-driven  quantum  system  near  or  adsorbed  on  a  surface, 
as  a  first  step  we  calculate  the  resonance  fluorescence  spectrum  of  a  two- 
level  atom  near  or  adsorbed  on  a  metal  surface.  In  the  surface-free  case, 
l.e..  In  the  absence  of  a  solid  surface  or  when  the  atom  Is  very  far  from 
the  surface,  resonance  excitation  of  a  two-level  atom  by  a  laser  field  has 

5 

been  extensively  investigated,  using  the  powerful  optical  Bloch  equations. 

There  has  also  been  Interest  In  surface-enhanced  spontaneous  emission  of 

6*9 

two-level  atoms  near  a  mirror.  In  this  paper,  we  provide  a  derivation 
of  the  surface-dressed  optical  Bloch  equations.  The  following  factors  are 
taken  Into  account: 

1)  Radiative  spontaneous  decay  and  stimulated  emission  produced  by  both 
the  driving  laser  field  and  the  reflected  electromagnetic  field. 

11)  Colllslonal  dephasing  of  the  two-level  atom  produced  by  for¬ 
eign  atoms  In  the  gas  medium,  which  constitutes  the  half 

region  of  the  Interface. 

Ill)  The  oscillatory  behavior  of  the  lifetime  of  the  adatom  due 
to  Interference  between  the  driving  field  and  the  reflected 


1 v)  Surface- Induced  dephasing  of  the  adatomic  transition  dipole, 
v)  The  Interaction  between  the  adatom  and  surface  plasmons. 
vl)  The  random  phase  fluctuation  of  the  laser  field  and  Its 
Influence  on  the  spectrum. 

In  Section  II  we  derive  the  surface-dressed  optical  Bloch  equations.  In 
Section  III  we  solve  the  equations  In  the  weak-fleld  limit  to  obtain  the 
population  Inversion  and  the  resonance  fluorescence  spectrum.  In  Section 
IV  we  present  results  for  these  quantities  and  the  surface- Induced 
phase-decay  constant  of  the  adatom.  Section  V  Is  the  Summary. 

II.  Surface-Dressed  Optical  Bloch  Equations  (SBE) 

Let  us  consider  a  laser-driven  two-level  atom  near  or  adsorbed  on  a 
metal  surface.  Although  the  atom  has  no  dipole  moment  In  Its  ground  state. 
It  can  have  a  significant  transition  dipole  connecting  the  ground  and 
excited  states.  The  reflected  field  provided  by  a  metallic  mirror  and 
surface  plasmon  resonance  can  Influence  the  dynamic  behavior  and  scat¬ 
tering  spectrum  of  the  adatom. 

To  treat  this  problem  we  utilize  a  self-consistent  approach.  This 
first  Involves  the  determination  of  the  Induced  transition  dipole  by 
quantum  mechanics,  l.e.,  the  optical  Bloch  equations.  The  transition 
dipole  Is  then  used  In  Maxwell's  equations  to  determine  the  reflected 
field,  and  the  reflected  field  In  turn  Is  used  In  the  quantum  mechanical 
procedure  to  find  the  transition  dipole.  The  resulting  form  of  the 
reflected  field  ER(t)  Is 


ER(t)  ■  vl2c12(t)f(d)  +  w21o21(t)f*(d). 


(1) 
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where  d  Is  the  distance  between  the  adatom  and  surface,  is  the 
electric-dipole  transition  matrix  element  between  the  states  |i>  and 
|j>.  £  | i ><j |  is  the  adatomic  transition  operator,  and  f(d)  is  a 

distance-dependent  function  which  will  be  discussed  in  more  detail 
later.  The  SBE  in  operator  form  are  then  given  by  (assuming  the 
equilibrium  condition  W  *  -1) 

£-02l(t)  *  (1«2rY2)°21(t)  +  *T**«(t)  *  [f(d)u12o12(t) 

♦  f*(d)u21o21(t)]> 


(2) 


*  -y[l+W(t)]  +  21^-  021(E(t)  ♦  Cf(d)p12a12(t) 

+  f*(d)ii21o21(t)]>  -2i  ol2(E(t)  +  [f(d)uj2o12(t) 


(3) 


♦  f*(d)u21o21(t)]} 


at 


012(t) 


(4) 


where  u>21  Is  the  adatomic  transition  frequency,  E(t)  is  the  random  phase- 
fluctuating  electric  field  strength  operator,  Yj  Is  the  population  relaxation 
constant,  y2  Is  the  phase  relaxation  constant,  and  fi(t)  s  ®22(t)  -  on(t)  Is  the 
population  Inversion  of  the  adatom.  Collisions  of  the  adatom  with  other 
atoms  In  the  gas  medium  are  accounted  for  in  the  Impact  limit  by  assuming 
Yj  f  2r2. 

The  quantum  equations  of  motion  for  W,  o12  and  o2l  are  most  useful  In 
the  rotating-wave  approximation  (RWA).  By  defining 


(S) 
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S21(t)  s 

where  Is  the  laser  frequency,  we  can  rewrite  Eq.  (2.2)  within  the 
MM  as 
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Q  (t)  s  nexp[+i$(t)]  Is  the  time-dependent  Rabl  frequency,  where  fi  s 
2 

Jrlu2ll E0^t)  and  E(t)  is  written  as  an  expectation  value  In  a 
coherent  state  of  the  laser  field  In  terms  of  the  phase  factor  *(t) 
as  E(t)  ■  EQ(t)exp[-1u^t  +  i$(t)]  +  c.c. ;  t  =  «2i"“L  ***•  detuning; 

the  surface- Induced  phase-decay  constant  Is  Im(f)|w12|2;  and 

^2  *  *2  +  V 

To  calculate  y*,  we  must  know  the  reflected  field  [Eq.  (1)]  or 
the  complex  function  f(d),  which  can  be  determined  by  the  Sommerfeld- 
Hertz  vector  procedure.2  For  the  case  of  a  laser-excited  atom  emitting 
near  a  metallic  half  space,  the  Imaginary  part  of  f(d)  has  the  following 
form: 
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ll 


(l+ej)^2 


2  “p 

b)  *  ■■*■■■ 

l+e, 


u462 

x  {— -  [(u2-u>2)2  +  ~^y  ~  e?u4][n  slnO  -  cosD] 
el  p  «  1  D* 

2  1 

+  2o>6u  [n  cosD  +  — k  sinD]}, 
p  D* 

where  c  ■  1  and  n  =  4c  +  i  for  the  case  of  the  induced  dipole 

DJ  D  , 

oriented  parallel  to  the  surface,  c  ■  2  and  n  =  -*  for  the  perpen- 

DJ 

dicular  case,  is  the  dielectric  constant  for  the  gas  medium,  kj  * 
wvTj/c  is  the  wave  number,  u  is  the  emission  frequency,  c  is  the  speed 
of  light,  and  D  ■  2kjd  is  the  reduced  distance  which  is  dimensionless. 
The  Drude  model10  has  been  used  to  determine  the  dielectric  constant 
of  the  metal  medium,  i.e.. 


(8) 


2 

e  (u)  «  i - -  ,  (9) 

u(ui  +  i  6) 

where  6  is  the  inverse  of  the  relaxation  time  and  up  Is  the  plasma 
frequency.  We  note  that  Im(f)  has  a  resonance  at  w._  ■  u>_(l  +  e,)"1^2, 

r  1 

which  Is  the  surface  pi asmon  frequency  of  the  Interface,  so  that  the 

reflected  field  given  by  Eq.  (1)  Involves  resonance  energy  transfer 

between  the  adatom  and  the  excited  surface  plasmon. 

Eq.  (8)  leads  to  a  surface-induced  phase-decay  constant  (in  the 

4u>3.  2 

unit  of  Einstein's  spontaneous  decay  constant  A  =  | u2i |  )  of  the 

following  form: 


II 
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y.  •  ca0[a,(n  slnO  -  -£2|2)  +  a?(n  cosD  +  £^)], 
*  u  1  D‘  i  Dz 


where 


a0  *1 


(1-e‘r  +  (a/u,Ap)‘ 


#1  -  (2-e2)2  +  (2a/^p)2  -  84 


a2  -  (28)  a/u.Ap,  (13] 

8  5  “L^p*  a  5  6/6.  C  and  n  are  the  same  as  In  Eq.  (8)  for  the 

parallel  and  perpendicular  cases,  and  c1  B  1  has  been  assumed  so 

that  w._  «  w/#  .  The  above  expression  for  y.  can  be  further 
P  t 

simplified  by  assuming  perfect  reflection,  i.e. ,  Rx  «  Rn  »  -1,  where  RU 
is  the  reflection  coefficient  for  an  Incident  ray  polarized  parallel  to 
the  plane  of  incidence  (p-polarized)  and  R1  is  the  reflection  coefficient 
for  perpendicular  polarization  (a-polarlzed).  Then 

\  c(n  sinD  -  (14) 


III.  Resonance  Fluroescence  Spectrum 

In  this  section  we  consider  the  weak  field  or  large  detuning 

limit,  where  W(t)  s  «W(t)»  *  -1.  We  also  use  the  phase-diffusion 
11  12 

model  ’  for  the  laser  field  In  order  to  Include  the  effects  of  the 
laser  bandwidth  yl:  • 


«n*(tl)n+(t2)»  »  nV^W. 


(15) 
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The  double  bracket  signifies  two  averages:  one  Is  over  the  stochastic 
ensemble,  and  the  other  Is  a  standard  quantum  mechanical  average.  Eq. 
(15)  applies  to  a  partial ly-coherent  laser  field,  which  Is  the  realistic 
situation.  Given  that  the  atom  near  the  metal  remains  more  or  less  In 
In  Its  ground  state,  we  can  linearize  Eq.  (7)  to  obtain 


$21^)  “  “(*2  +  “  *A)S2i(t)  0  (t). 


(16) 


To  determine  the  population  Inversion,  let  us  take  the  laser 
field  to  be  turned  on  at  t  ■  0  with  a  constant  amplitude  for 
t  ^  0.  We  then  obtain  the  analytical  form 


W(t)  «  -1  +  n  { 


2,  V*L 


n[(VrLra 


yi‘Vyl 

—  — 7 — r 
Yl^Yl“Y2“YL^  +a * 


t(VYL^Yl“VYl.)  +  A23««Ut)  ♦a(2^2-2YL-Y1)s1n(at) 
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x  e’^2  + 


(17) 


To  obtain  the  power  spectrum  of  the  scattered  light,  we  must  find 

A  A 

the  dipole-dipole  correlation  function  «S21(t2)S12(tj)».  For 
tl*t2>>  y11,y21  th1s  *****  the  fonn  (t2  I  V 


<<s2i(t2)Si2(ti)>> 


n  l 

^  <y,,W* 


.  t(vy i><yv<»HiKy\>-u3  ,(t  -uh^-v 


(continued) 
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(continued) 


e-YL(t2-t1) 


+  T  (?2*YL+U)(?2-YL-Uj 


(18) 


The  corresponding  power  spectrum  for  the  emitted  frequency  u  Is 
given  by  the  Fourier  transform  of  the  above  equation: 


s(«)  « -V 


(>2+Yl)2+A2  (Y2*YL)t+A' 


1  Y|_  (^2“  yl+a^  ^  +  Zy^w-u^) 

TT“2  t— - 


T77 — nr 

Y^  +  (w-u>^ ) 


f2(f|-Tf*42)  +  2yl«(u-»21),  2(f2»YL)f? 

- - 2 - - J  +  77  :  .77 ]  v  ; 

Y2  +  (a)-ai2j! 


Y9  +  («,,?  Yj[?2  +  (-W^fi 


The  fluorescence  spectrum  of  Eq.  (19)  generally  exhibits  two  peaks 
for  d  surface-dressed  two-level  atom  In  a  weak,  partial ly-coherent 
laser  field.  One  peak  Is  centered  at  «  •  which  corresponds  to 
elastic  or  Rayleigh  scattering,  and  the  other  near  u  ■  «2j  corresponds 
to  Inelastic  scattering  or  fluorescence.  In  the  surface-free  case  for 
a  colllslonltss  atom  In  a  ful ly-coherent  laser  field  (yj  *  2y2  »  A, 
yl  »  0),  the  above  spectrum  becomes  a  single  peak. 


S(«) 


*P2/2 

(A/2)Sa2 


6(u-w^), 


(20) 
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which  Is  the  standard  result  of  Heitler.  In  general,  however,  the 
adatom  suffers  collisions  with  foreign  atoms  In  the  gas  medium,  so 
that  y  f  2y2>  and  there  are  nonzero  values  for  the  surface-induced 
phase-decay  constant  y t  and  the  laser  bandwidth  yL«  This  Is  then 
responsible  for  the  second  Lorentzlan  peak  occurring  near  the  transl 
tlon  frequency  u>21. 

IV.  Results 

Fig.  1  shows  the  time  evolution  of  the  population  inversion  W(t) 
[Eq.  (17)]  for  an  atom  near  a  metal  surface  for  two  different  adatom 
surface  plasmon  resonance  conditions,  8  =  w.  /w._  ■  0.5  and  b  =  0.9. 

L 

-3 

We  have  used  the  value  of  3.5  x  10  for  the  ratio  6/w._,  which 
corresponds  to  a  silver  surface.  Each  curve  approaches  a  stationary 
state  value,  assuming  that  the  duration  of  the  laser  pulse  is  much 
longer  than  the  relaxation  time  of  the  adatom  (characterized  by  y^1 
and  ?2*).  The  reason  the  values  of  W  for  the  off-resonance  case 
(8  *  0.5)  are  less  than  those  for  the  near-resonance  case  (e  *  0.9) 
is  that  the  levels  of  the  adatom  in  the  latter  case  are  strongly 
broadened  by  the  resonant  transfer  of  energy  between  the  adatom 
and  surface  plasmon,  so  that  the  population  of  the  upper  level 
increases.  The  oscillations  increase  and  the  damping  due  to  the 
reflected  field  disappears  as  the  adatom  moves  further  from  the 
surface  (going  from  curve  1  to  3).  The  strong  oscillations  in  the 
nearly  surface-free  case  (curve  3)  are  caused  by  beating  between  the 
driving  laser  field  and  the  "adatomic  oscillator". 

Fig.  2  displays  the  resonance  fluorescence  spectrum  [Eq.  (19)], 
with  the  Rayleigh  and  fluorescence  peaks.  The  former  peak  in  our 
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model  Is  just  determined  by  the  laser  bandwidth  y^.  The  latter  peak  is 
strongly  dependent  on  the  adatom-surface  distance;  as  the  distance 
decreases,  the  upper  adatomic  level  broadens  and  hence  the  fluroescence 
intensity  increases. 

In  Fig.  3  we  plot  W(t)  for  the  case  of  a  perfectly  reflective 
metal  surface,  R"^*  R11  *  -1  [see  Eq.  (14)],  for  several  values  of  the 
laser  bandwidth  and  adatom-surface  distance.  The  oscillations  increase 
as  the  adatom  moves  further  from  the  surface  and  decrease  as  the  laser 
bandwidth  increases  [going  from  (a)  to  (c)].  A  larger  bandwidth  leads 
to  more  damping  of  the  adatomic  oscillator  [see  Eq.  (17)]  and  hence 
an  increase  in  the  population  inversion.  It  is  interesting  to  note 
that  W  is  more  sensitive  to  the  adatom-surface  distance  for  the  case 
of  the  transition  dipole  oriented  parallel  to  the  surface  than  for  the 
perpendicular  case.  We  see  that  the  positions  of  the  dashed  curves  2 
and  3  are  interchanged,  indicating  that  the  power  spectrum  oscillates  as 
a  function  of  D.  As  seen  by  the  form  of  Eq.  (8),  we  would  also  see  such 
oscillations  for  the  population  inversion  and  power  spectrum  in  Figs.  1 
and  2  if  we  were  to  consider  a  wider  range  of  values  of  D. 

Fig.  4  displays  the  surface-induced  phase-decay  constant  y.  as  a 
function  of  the  adatom-surface  distance  for  a  perfectly  reflective 
surface.  The  oscillations  are  due  to  the  Interference  between  the 
reflected  field  and  the  driving  field.  This  gives  rise  to  oscillations 
in  the  lifetime  of  the  adatomic  oscillator,  while  the  small  negative  value 
of  yt  in  some  regions  of  0  tends  to  Increase  the  lifetime. 

Fig.  5  shows  the  power  spectrum  for  several  values  of  the  laser  band¬ 
width  yl-  When  yL  is  very  small  compared  to  the  phase-decay  constant  yj 
[Fig.  5(a)], 
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the  power  Is  concentrated  largely  In  the  Rayleigh  peak.  Here  the  population 
of  the  upper  adatomlc  level,  and  hence  the  power  of  the  fluorescence,  is 
determined  mainly  by  the  broadening  of  that  level  due  to  the  surface,  so 
that  the  power  Increases  as  0  decreases.  The  situation  Is  changed  when 
yl  Is  larger  than  y2»  In  which  case  the  power  of  the  Rayleigh  component 
transfers  to  the  fluorescence  peak  [Fig.  5(b)  and  (c)].  Although  a  decrease 
In  D  leads  to  a  broadening  of  the  adatomlc  level,  now  the  total  power  in  the 
fluorescence  peak  Is  approximately  unchanged,  whereby  a  curve  for  smaller  D 
has  a  larger  width  but  a  smaller  peak  height.  A  more  detailed  discussion  on 
the  relation  between  fluorescence  and  laser  bandwidth  in  general  will  be 
published  later. 

It  is  interesting  to  note  that  Weitz  et  al_  recently  published  mea¬ 
surements  and  a  theoretical  model  on  a  hierarchy  of  enhancement  ratios  for 

normal  Raman  scattering,  resonance  Raman  scattering  and  fluorescence  from 

14 

molecules  adsorbed  on  a  silver  surface;  The  results  from  our  two- level 
adatom  model  for  surface-enhanced  fluorescence  are  in  good  agreement  with 
their  measured  enhancement  ratio  of  10” 1  to  10  for  fluorescence. 

V.  Summary 

A  simple  model  has  been  developed  to  describe  the  dynamic  behavior 
of  a  two-level  adatom  at  a  metal  surface.  Surface-dressed  optical  Bloch 
equations  (SBE),  which  account  for  surface- reflected  photons  and  adatom- 
surface  plasmon  coupling,  have  been  derived.  Collisions  with  foreign  gas 
atoms  and  effects  of  the  laser  bandwidth  have  been  Included.  The  SBE  have 
been  used  to  evaluate  the  power  spectrum  of  the  scattered  light,  the  popula¬ 
tion  Inversion  and  the  surface- induced  phase-decay  constant  as  functions  of 
the  adatom- surface  distance  and  adatom-surface  plasmon  resonance  parameter 
B.  The  displayed  properties  of  the  laser-adatom- surface  system  should  be 
useful  In  understanding  basic  processes  occurring  In  surface  resonance 
spectroscopy. 
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Figure  Captions 

1.  Time  evolution  of  the  population  Inversion  In  the  weak  field  or 
large  detuning  limit,  with  (y^t  a.  a)  •  (0.3.  0.05,  5) 

In  the  unit  A  for  0  »  0.9  (solid  curves)  and  0  ■  0.5  (dashed 
curves).  Curves  1,  2  and  3  correspond  to  the  values  of  the 
reduced  distance  0  ■  2kjd  of  1,  3  and  5,  respectively.  The 
horizontal  axis  corresponds  to  time  In  the  unit  A’*.  The  quantity 
[W(t)  +  l]/o2  Is  shown  along  the  vertical  axis. 

2.  Fluorescence  power  spectrum  of  scattered  light  with  (yl.  Q,  a)  ■  (0.3, 
0.05,  5)  in  the  unit  A  and  0  ■  0.9  for  the  Induced  transition. dipole 
oriented  perpendicular  to  the  surface.  The  units  along  the  vertical 
axis  are  arbitrary.  Curves  1,  2  and  3  correspond  to  0  ■  1,  3  and  5, 

•  respectively. 

3.  Time  evolution  of  the  population  Inversion  W(t)  of  an  adatom  at  a 
perfectly  reflective  metal  surface,  with  Q  *  0.05,  a  »  5,  (a)  yl  * 

0.3,  (b)  yl  *  1,  and  (c)  yl  *  3  (In  the  unit  A).  The  solid  curves 
are  for  the  Induced  transition  dipole  oriented  perpendicular  to 
the  surface,  while  the  dashed  curves  are  for  the  parallel  case,  and 
the  labels  1,  2  and  3  correspond  to  0  ■  2,  4  and  6,  respectively. 

4.  Surface-Induced  phase-decay  constant  y4  (In  the  unit  A)  as  a  function 
of  the  reduced  adatom-surface  distance  D  ■  2kjd  for  a  perfectly 
reflective  metal  surface.  The  solid  curve  Is  for  the  Induced 
transition  dipole  oriented  perpendicular  to  the  surface,  while 

the  dashed  curve  Is  for  the  parallel  case. 

5.  Power  spectrum  of  scattered  light  for  a  perfectly  reflective  metal 
surface,  with  n  ■  0.05,  a  ■  5,  (a)  yl  •  0.3,  (b)  yl  ■  1  and  (c)  yl  • 

3  (In  the  unit  A).  The  Induced  transition  dipole  Is  oriented  perpendi¬ 
cular  to  the  surface.  The  units  along  the  vertical  axis  are  arbitrary. 
Curves  1,  2  and  3  correspond  to  0  •  2,  4  end  6,  respectively. 
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